home *** CD-ROM | disk | FTP | other *** search
/ Ham Radio 2000 #2 / Ham Radio 2000 - Volume 2.iso / HAMV2 / ANTENNA / YAGIU112 / MPR_HACK.C < prev    next >
C/C++ Source or Header  |  1995-08-07  |  324b  |  19 lines
  1. #include "yagi.h"
  2.  
  3. void mprove(double **a,double **alud,int n,int *indx,double *b,double *x)
  4. {
  5.     int j,i;
  6.     double sdp;
  7.     double *r;
  8.  
  9.     r=dvector(1,n);
  10.     for (i=1;i<=n;i++) {
  11.         sdp = -b[i];
  12.         for (j=1;j<=n;j++) sdp += a[i][j]*x[j];
  13.         r[i]=sdp;
  14.     }
  15.     lubksb(alud,n,indx,r);
  16.     for (i=1;i<=n;i++) x[i] -= r[i];
  17.     free_dvector(r,1,n);
  18. }
  19.